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The generalized t-J model conserving the number of double occupancies is constructed from the 
Hubbard model at and in the vicinity of half-filling at strong coupling. The construction is realized 
by a self-similar continuous unitary transformation. The flow equation is closed by a truncation 
scheme based on the spatial range of processes. We analyze the conditions under which the t-J 
model can be set up and we find that it can only be defined for sufficiently large interaction. There, 
the parameters of the effective model are determined. 
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I. INTRODUCTION 

The Hubbard model serves as a prototype model for 
the description of strongly correlated electron systems. 
It consists of electrons with spin moving on a lattice. 
They repel each other on-site which leads to a complex 
interplay of magnetic and charge degrees of freedom. At 
half-filling, the model displays insulating behavior in the 
limit of strong interaction. In this limit, the model can 
be mapped onto a pure spin Heisenberg mode l 1 ! 2 ! 3 ^ 5 ! 6 . 
On the other side, vanishing interaction makes it a metal. 
In-between an intricate metal-insulator transition (MIT) 
takes place which is well understood only in extreme 
cases. In one dimension, the Bethe ansatz tells us that ar- 
bitrarily small interactions render the model insulating'. 
In infinite dimensions, the dynamical mean-field the- 
ory allows to make statements about the value of the 
critical interaction 8-910 if long-range magnetic order is 
suppressed. It is argued that this suppression can be 
achieved physically by introducing frustration. 

It is the aim of the present paper to provide a system- 
atic and controlled non-perturbative derivation of the t-J 
model from the Hubbard model. Most importantly, we 
will discuss under which circumstances the reduction of 
the Hubbard model to a t-J model is justified. The t-J 
model is the effective model describing fermions which 
interact and hop without creating or annihilating dou- 
ble occupancies (DOs). This means that never two elec- 
trons occur on the same site. Thus the reduction to a 
t-J model corresponds to the elimination of charge fluc- 
tuations. The parameters of the effective model will be 
computed quantitatively. Our results reach in two ways 
beyond what has been done beforei*2*ii^&. First, we 
provide non-perturbative results. Second, we discuss the 
matrix elements of hole motion and interaction in a sys- 
tematic way beyond the zeroth order. Both ingredients 
help us to discuss the breakdown of the reduction to a 
t-J model. 

The reference ensemble in our approach is the magnet- 
ically completely disordered half-filled model (see Sect. 
III). This constitutes the vacuum in our calculations. In 
this sense, we start from the Hubbard model at half- 
filling. But our approach naturally keeps also track of 
the dynamics of holes. This is unavoidable since we have 



to know how virtual intermediate states evolve and they 
contain generically charge excitations. The parts of the 
effective model which describe the dynamics of holes ap- 
ply also to holes inserted externally by doping. In this 
sense, the effective model describes also the physics of a 
vicinity of half-filling. But it is beyond the scope of the 
present paper to discuss the doping dependence of the 
couplings of the effective model since this would corre- 
spond to a change of the reference ensemble. 

In Fig. ^ a schematic phase diagram is shown which 





1 T 

/ 


1 1 


PI 


/ 


PM 




AFI 





1 2 



w/u 

FIG. 1: Schematic phase diagram of the Hubbard model at 
half-filling in higher dimensions d > 2. The tick labels at the 
x-axis are indicative only. The bare band width is denoted 
by W, the Hubbard repulsion by U, the temperature by T; 
PI stands for paramagnetic insulator, PM for paramagnetic 
metal and AFI for antiferromagnetic insulator. The dashed 
line indicates the transition or the crossover from insulator to 
metal. 

summarizes the currently assumed picture. This picture 
is motivated from various mean-field computations, in 
particular in infinite dimensions^. For large repulsion 
U (left side in the figure), the electrons are localized so 
that the system is insulating. The spins interact via the 
magnetic exchange J of the order of At 2 /U which sets 
the energy scale for the transition between a long-range 
magnetically ordered phase (AFI) and a disordered phase 
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(PI). For decreasing U, the ground state remains insulat- 
ing unless the magnetic order is fully suppressed. If the 
spin state becomes sufficiently disordered, for instance 
at higher temperatures, the system becomes conducting 
(PM). The change from paramagnetic insulator to para- 
magnetic metal is a first order transition in infinite di- 
mensions and for not too large temperatures^. In gen- 
eral, it is to be expected that the system undergoes a 
crossover. In two dimensions in particular, the case in- 
vestigated here, no long-range antiferromagnetism occurs 
at finite temperaturesii so that the AFI collapses and 
only the disordered phases occur. 

Once the interaction is such that the system is a para- 
magnetic metal, charge excitations are possible at zero 
energy. Hence, it cannot be expected that the elimina- 
tion of charge fluctuations, virtual double occupancies, 
is possible. This means that the reduction of the Hub- 
bard model to a t-J model can only be defined as long 
as the paramagnetic phase is insulating. The sole oc- 
curence of an insulating phase with long-range order is 
not sufficient. The insulating behavior in the long-range 
ordered AFI phase for interactions where the paramag- 
netic phase is metallic is due to the concomitant breaking 
of the translation symmetry. Hence it resembles rather a 
band insulator than a genuine Mott-Hubbard insulator—. 
The t-J model, however, describes low-lying degrees of 
freedom, namely the degrees of freedom of the spins and 
of doped charge carriers, without the precondition of a 
certain long-range order. The derivation of the t-J model 
requires to eliminate the virtual charge fluctuations while 
leaving the state of the spins and of the doped charges 
unspecified. One may consider the unspecified low-lying 
degrees of freedom to be at infinite temperature in a com- 
pletely disordered mixture. The phase diagram in Fig. 2] 
tells us that there must be a certain U of the order of the 
band width W below which the Hubbard model cannot 
be represented by a t-J model. 

The method of continuous unitary transformations 
(CUTs) has been proposed by Wegner in the context of 
condensed matter physics^. Simultaneously, similar ap- 
proaches were tested by Glazek and Wilson in high en- 
ergy physics^^. The proper choice of an infinitesimal 
generator brings the many-body problem under study 
into a more tractable form by decoupling different sec- 
tors of the Hamiltonian. In many cases it is possible to 
advance to an effective model that is tractable by some 
other means while retaining the complex physics of the 
original problem 6 i 16 i 17 i 18 i 19 i 20 . We will apply a CUT to 
the Hubbard model in order to eliminate operators that 
couple states of different double occupancy. This has 
been done before by Stein& in a perturbative fashion. In 
contrast to Stein's perturbative treatment we will adopt a 
self-similar truncation scheme for the operators emerging 
during the transformation. This self-similar CUT retains 
terms of a certain structure. 

The paper is organized as follows. The Hubbard model 
is introduced in Sect. [H] Sect. IIIII explains the method 
of self-similar CUTs. A transparent example is given to 



illustrate technical details and characteristic properties 
of the CUT applied to the Hubbard model. The results 
for the range of validity of the mapping and for the ef- 
fective parameters are presented in Sect. lIVl Evidence is 
provided that the mapping cannot be defined for all ra- 
tios W/U . The physical implications of our findings are 
discussed in Sect. [3 A brief summary is given finally in 
Sect-Ell 



II. HUBBARD MODEL AND OPERATORS 

Our starting point is the Hubbard model on a two- 
dimensional square lattice in the vicinity of half-filling 
with nearest-neighbor hopping. Since the calculation 
covers also the motion and interaction of charge carri- 
ers it is not limited to strictly half-filling but covers also 
small hole concentrations. 

The Hamilton operator is split into a kinetic part H t 
and an interaction Hjj 

H = Hu+H t (la) 
H v = Uj2(ni,i-l/2)(n u -l/2) (lb) 

i 

= (U/2)(b-N/2) 

H * = 1 E ( c L^+ h - c -). ( lc ) 

where N is the number of sites, c icr> c i are creation and 
annihilation operators of an electron on site i with spin 
a E {T>4} and n iy(7 — c\ rr c ic is the number operator. 
Hopping is only possible between nearest neighbor (NN) 
sites as indicated by in the sum in H t . The bare 

hopping coefficient is given by t so that in two dimensions 
the band width is W = St. It turns out to be convenient 
to denote all results in terms of the band width W. 

We also define the operator D :— Y^i[ n iA n i,l + (1 — 
rii^)(l — rii t i)] that counts the number of double occu- 
pancies (DOs) of particles and holes. Hjj describes the 
repulsion of electrons of opposite spin on the same site. 
For large U the density of states splits into an upper and 
a lower Hubbard band separated by an energy of order U. 
In the half-filled case, in the limit of infinite U the elec- 
trons are fixed on their lattice site by the constraint of 
having no DO. Reducing U from this limit the electrons 
are gradually allowed to move producing DOs intermedi- 
ately. But the physics remains still rather local. In the 
following, we want to map the Hubbard model for finite 
U in a systematic way onto an effective model that con- 
serves the number of DOs. The guiding idea will be that 
the physics can be incorporated into an effective model 
that contains operators that describe rather local pro- 
cesses. For too small values of U the vanishing of the 
charge gap will make the mapping impossible. 
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III. CONTINUOUS UNITARY 
TRANSFORMATION 

In this section, the technical details of our calculations 
are presented. The method of continuous unitary trans- 
formations (CUT) introduces a continuous auxiliary vari- 
able £ and transforms the Hamiltonian according to the 
flow equation^ 



df 



H(£) = [r,(£),H(£)] 



(2) 



with an antihcrmitian generator n(£). The operators 
H (£) for different values of £ are unitarily equivalent. The 
transformation starts with the initial condition H(£ = 
0) = H and ends with an effective Hamiltonian at £ — > oo. 
We want the effective Hamiltonian to conserve the num- 
ber of DOs and therefore use the generator— 



r,(£) = [D,H(£)] 



(3) 



Except for an ov erall facto»2i, the above 77 coincides with 
r\ in Refs. Il7l22t where it was defined using the sign of 
the change in the number of quasiparticles ensuring that 
the block-band structure of the Hamiltonian is conserved. 
This means that during the flow no operators are gener- 
ated that change the number of DOs by a value other 
than 0, -2 and 2. 

To see how the CUT works, it is helpful to classify the 
operators in the kinetic part of the hamiltonian according 
to their effect on the number of DOs 



H t = T +T +2 + T. 



Tn = t 



(4a) 



(4b) 



T+2 = t+2 E [ n l ,<T C U C jA 1 - n j.a) 



>c I-q, „ (1 - n-) (4c) 



+ C 1 - n j,M,<? c i,w n i,°\ > ( 4d ) 



where a = —o~. The projection operators ensure that T n 
changes the number of DOs by n. 

Operators will be denoted in a standardized normal- 
ordered real space representation (see Table QJ. The pre- 
factors of such a product of local operators will be called 
coefficient of this term in the Hamiltonian. Eq. will 
yield differential equations for these coefficients. The co- 
efficients to, t+2 and t^2 coincide and are equal to the 
bare hopping t in the beginning of the transformation. 
But they will develop differently during the flow. 



The operator T n yields a contribution [D, T n (£)] 
nT n {£) to n leading to 



= -^T n {t) 



(5) 



and thus to a suppression of any term not conserving the 
number of DOs. 

In general the prescription (0) will not produce closed 
equations. One has to find a way to deal with the pro- 
liferating number of terms. Stein-i set up a perturbative 
CUT for the Hubbard model with the same generator J3J). 
He ordered the terms produced on the right hand side of 
(J2J) according to the power oit/U with which they were 
generated. The terms were kept in a form consisting of 
long products of To, T+2 and T_2- 

Contrary to this perturbative strategy, we will not keep 
the commutators of To, T+ 2 and T_2 in the unevalu- 
ated form. We compute the right hand side of the flow 
equation @ explicitly. This will produce terms already 
present in the Hamiltonian and new terms. According to 
systematic rules the new terms will be kept or discarded. 
Finally, a closed set of terms is reached. The differential 
equations representing the flow equations of the CUT are 
obtained by comparing the coefficients of all the retained 
terms. Since both n and H are linear in the coefficients 
the right hand sides of the differential equations are bi- 
linear in the coefficients in close analogy to conventional 
renormalization group equations. 

To find a scheme to truncate the generation of new 
terms we make the following considerations. In the limit 
of large U the motion of the electrons at half-filling is 
suppressed by the energy cost of DOs. Therefore the op- 
erators that become important first when going to lower 
U are those that describe local processes. To use this for 
the truncation scheme we have to define a measure for 
the locality of a term. An undispensable prerequisite is 
a physically meaningful and in particular unique way to 
denote the operators. Hence we define a sort of normal- 
ordering of the operator products^. We will normal- 
order all terms with respect to the half-filled paramag- 
netic reference ensemble represented by the statistical op- 
erator 



Po 



n(i/2)(IT)«<T| + U)«UI) 



(6) 



where we use a basis of local states at each site i. There 
are four local states: the empty site |0), two singly oc- 
cupied states with spin up und down, | |) and | j), and 
the state with two electrons, | fj.). The reference ensem- 
ble © contains only states with one electron per site 
since it is half-filled. The up-spin and down-spin states 
are weighted equally to represent a completely disordered 
paramagnetic reference ensemble. Any deviation from 
the reference ensemble represents an excitation. 

A product of n local operators 0\ ■ ... ■ O n is normal- 
ordered if any expectation value 



(Oi) lei := Tr(O lPQ ) 

= (1/2)(,(T|C,|T> 1 



i(i\Oi\i)i) 



(7a) 
(7b) 
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with the reference ensemble vanishes. In Table[I]all local 
operators are given in their normal-ordered form. Apart 
from the unity matrix 1 all local operators in Tab. [I] are 
chosen such that the expectation value Q with the refer- 
ence ensemble yields zero. The operator n = n* + n, — 1 
counts the number of electrons relative to half-filling. 
The Hubbard operator 2n^n^ — n counts the number of 
DOs. The empty state |0) and the doubly occupied state 
| U) yield (0|2n T n i - n\0) = (U - n\ U) = 1 

whereas expectation values with 1 1) and 1 1) vanish. For 
the spin operator a z the sum {7J is zero as well. All other 
operators in Table [I] are non-diagonal in the local basis 
and thus normal-ordered with respect to the reference 
ensemble. 



G Z = 


1 


n = + — 
2n^n^ — n 


1 


g+ = c\ Cl 

a' = C |c T 


C | C T 
C T C I 


(1 
(1 


-n | )c T 


(1 - n T )c x 
(l-n T )cJ 




n A 


n i c T 



TABLE I: Local operators. The operators in the first two 
lines have an even number of fermionic operators, those in 
the last two lines an odd number of fermionic operators; the 
usual spin operators are S z — g z /2, S — a . 

The local operators of Tab.[I]have to be labelled with a 
site-index when they are used to write down the Hamil- 
tonian. Consequently, also products of local normal- 
ordered operators belonging to different sites are normal- 
ordered. The expectation values of each separate factor 
with respect to the reference ensemble vanish. It is the 
prefactor of a product of local operators as defined in 
Tab. HI which we will call the coefficient of this term. 
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FIG. 2: The operator from Eq.[H]in the 2-dimensional square 
lattice. Filled sites are involved in the hopping process. The 
operators are shown at the sites on which they act. 



this process, which is compensated by the spin-flip on site 
i. This term is illustrated in Fig. |21 It contains operators 
that act on three different sites. They are shown as filled 
circles in Fig. El These sites will be called the cluster be- 
longing to this term. The maximal taxi cab distance d of 
two sites in this cluster is the extension of this term. The 
extension of the term © in Figgis d = 2. The extension 
of the cluster belonging to a term will be our measure of 
its locality. For a truncation scheme, we define a certain 
set of clusters, i.e. geometric arrangements of sites which 
are translated over the whole lattice. Those operators 
generated in the flow equation @ are kept which can be 
defined completely on one of the clusters. If they affect 
more sites or sites at a larger distance the operators are 
discarded. 

We emphasize that the use of a set of clusters to define 
the truncation does not turn our approach into a finite 
cluster calculation. Rather it is a self-similar, renormal- 
izing calculation in the thermodynamic limit where the 
decision which terms are kept is based on their locality. 
By translational invariance, the operators appear in the 
whole lattice as is denoted by the J2i m Eq. ©• F° r large 
enough clusters, the perturbative results are reproduced 
by the self-similar calculation. For simplicity, the differ- 
ential equation is set up for the coefficient of a single rep- 
resentative of a certain kind of term taking advantage of 
the symmetries of the problem, namely translational in- 
variance, the point-group symmetries, spin rotation sym- 
metry, and hermiticity. 

The complete prescription for the self similar CUT 
calculation reads: (i) Define a certain set of clusters. 
Normal-ordered operators on these clusters are kept dur- 
ing the calculation. Additional conditions may be im- 
posed. With this truncation rule the flow equation @ 
closes, (ii) Calculate [n(£), H(£)] to obtain the flow equa- 
tion J2J) for the coefficients of the terms in H. Terms 
not present in the original Hamiltonian start at £ = 
with initial coefficients zero, (iii) The numerical integra- 
tion of the flow equation yields the effective model for 
£ — > oo; it conserves the number of DOs. This step re- 
quires that the flow equation converges for £ — > oo. If it 
converges the effective model is constructed successfully. 
Generally, non-convergence can be due to an insufficient 
approximation or due to the breakdown of the mapping 
of the Hubbard model to the t-J model. 

Our procedure for a self similar CUT calculation will 
be exemplified in the following subsection. 



On the basis of the normal-ordering and of the stan- 
dardized notation introduced above we define a measure 
for the locality of a term, namely its extension d. Con- 
sider the following exemplary normal-ordered term 

X! n i+e y ,1 C i+e y ,l c i+e !c ,1 n i+e x ,l ' a i '• ( 8 ) 
i 

This term represents the hopping of an electron from site 
i + e x to site i + e y . The spin of the electron is flipped in 



A. Example: The minimal and the NN truncation 

Here we work out a simple example and compare its 
results to perturbation theory. At first we build a mini- 
mal model for the NN coupling J. Starting with H t and 
Hjj the generator r] reads 

T}(£) - [D, H{£)] = [D, T +2 + T_ 2 ] = 2X+ 2 - 2Tl 2 (9) 
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To set up the flow equation we have to calculate 
d 



df: 



H{£) = WlH{i)\ = W),Hu(£) + H t (£)]. (10) 



The first contribution to the flow equation 

[r,(l), Hull)} = -2U(£)T +2 - 2U(£)T- 2 (11) 

produces the suppression of terms changing the number 
of DOs as intended. The second contribution reads 



[ V (£),H t (£)} = [2T +2 -2T_ 2 ,T +T +2 + T_ 2 ] 



(12a) 
(12b) 



+ 16t +2 t_ 2 S * ' S J 

(id) 

UiUj (12c) 

(id) 

+st +2 t_ 2 £ cyy u c jA (i2d) 

(id) 

+... (three site terms). 

In Eq. 112|) the terms containing operators on three dif- 
ferent sites are omitted. The term (I12at will renormalize 
the strength of the Hubbard interaction U. The next line 
(|12b(l generates the Heisenberg exchange 



i?NN — Ji(£) ^ Si • Sj 



(13) 



where C = arctanh(2J7o/s). The effective model is ob- 
tained in the limit £ — > oo. Since t +2 {oo) = it contains 
only #DO-conserving terms 



H cS = To + U{<x)\b + Ji(oo) S * ■ S J 



(16) 



(id) 



Strictly at half-filling, Eq. 1|16|) reduces to a Heisenberg 
model with effective coupling Ji je fi = Ji(oo) 



Ji,cff - ^\/l + lS(t/U) 2 -^ (17a) 



U 
At 2 



o(t 4 /u 3 ). 



(17b) 



For small t/U this reproduces the well-known second or- 
der result 4 — ^jj-. Note that the final results are 
given in t and U denoting the initial, bare values whereas, 
for clarity, we use to = t and Uq — U to denote the initial, 
bare values while solving the flow equation. 
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FIG. 3: Filled circles indicate the two clusters defining the 
truncation rule of the NN truncation. Only terms fitting on 
these clusters are kept. 



■d) 



with the initial condition Ji(0) = 0. To get a minimal 
model for the NN exchange we neglect the terms in (I12cf) 
and l|12d() . The exchange term i?NN will not produce a 
contribution to rj since it does not change the number 
of DOs. As a further simplification we neglect the terms 
that would arise from the commutator [rj, Ann]- Exploit- 
ing that operators related by hermitian conjugation have 
the same coefficient we know t +2 = f_ 2 . Thus the differ- 
ential equations for the minimal truncation read 



d e U{£) = 32t +2 (£) 2 

d t t (£) = 

d f t +2 (£) = -2U(£)t +2 (£) 

d t j x {£) = m +2 {£f , 



(14a) 
(14b) 
(14c) 
(14d) 



where we use dg as shorthand for 4p. Using the con- 
served quantity s = y/ 64t +2 (£) 2 + 4U (£) 2 we obtain for 
the minimal truncation 



U(£) = -tanh(s^ + C) 



(15a) 

t (£) = to(0) (15b) 

t +2 (l) = t_ 2 (£) = (s/8)yl - tanh 2 (s^ + C) (15c) 
Ji(£) = s [tanh(s£ + C) - 2U /s] /4 (15d) 



The result l|17(l coincides with the result one gets from 
diagonalization of a two-site cluster for the splitting be- 
tween the singlet and the triplet state and thus for J. 
This is purely coincidental since taking into account all 
operators on NN sites in the CUT leads to a different 
result. The truncation rule for this NN truncation is de- 
fined by the two clusters shown in Fig. |3| All terms are 
kept that contain only local operators on two neighboring 
sites. For the effective coupling we find, see Appendix lAl 



NN 
l.cff 



1 + 28 (t/Uf 



(18) 



As expected the result is different from the result for two 
sites since the exact diagonalization deals with a finite 
system only whereas the CUT is a renormalization pro- 
cedure on the thermodynamic lattice where the terms 
kept are chosen according to their locality. 

To compare our results to perturbation theory, we cal- 
culated also the 4th and 6th order in t/U for J\ in per- 
turbation theorjiiiSi. Fig. 0] compares the value of the 
Heisenberg coupling J\ from the minimal and from the 
NN truncation to the results from perturbation theory. 
All results are given relative to the second order result 
= 4t 2 /U. The perturbative results show strong di- 
vergencies even if one includes the 4th and the 6th order. 
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FIG. 4: Results for the Heisenberg coupling Ji from the mini- 
mal and the NN truncation compared to perturbation theory. 
All curves are plotted relative to j[ = 4t 2 /U. The numbers 
in the legend in front of Pade or DlogP indicate the degree of 
the polynomials in x — t 2 /U 2 in the numerator and denomi- 
nator of the Pade approximant. Recall W = 8t. 



neighboring sites is sufficient. Thus the NN truncation 
is defined by the cluster given in Fig. [SJi. The cluster in 
Fig-EK represents the two possible clusters of NN sites on 
the two-dimensional square lattice. In 4th order in the 
electron hopping t also further hopping processes con- 
tribute to the effective spin-couplings. The cluster of four 
sites on a plaquette of the square lattice and the linear 
clusters of three sites cover all relevant hopping processes 
in 4th order. They are depicted in Fig. [5p. The double 
plaquette and the linear chain of four sites depicted in 
Fig. are representatives for the clusters necessary for 
accounting fully of the 6th order. 
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The dashed and solid black lines show the solution for J\ 
from the minimal and NN truncation, respectively. They 
give meaningful finite values for J\ without divergent be- 
havior up to high values of W/U. Already the minimal 
model correctly reproduces the leading perturbative be- 
havior oc At 2 /U for small W/U. 

In addition, Fig. 0] shows Pade-approximants derived 
from the series of the perturbation theory. The contin- 
uous grey line depicts a Dlog-Padc- approximant. The 
Pade-approximants support that the CUT results for J\ 
lie in the reasonable range of values. But neither the 
CUT nor the Pade results determine J\ reliably for values 
of W/U > 0.8. Therefore, we improve the CUT calcula- 
tion by taking into account more terms on more complex 
site configurations leading to higher truncation schemes. 
The agreement of calculations with various numbers of 
operators will be a probe for the accuracy of the results. 



B. Higher truncation schemes 

Fig. shows the clusters defining the various trunca- 
tion schemes. The sites that belong to the clusters are de- 
picted as filled circles within the two-dimensional square 
lattice. Only one representative is shown for different 
clusters that are related by the point-group symmetry 
of the square lattice. We use the clusters that are neces- 
sary for a perturbative derivation of the spin Hamiltonian 
from the Hubbard model as a guideline for the choice 
of the clusters defining the truncation. In second order 
the NN Heisenberg exchange is the only effective spin- 
coupling. For its perturbative derivation a cluster of two 



FIG. 5: Various clusters defining the operators to be kept in 
the CUT. Only one of various clusters that are related by the 
point group symmetry of the two-dimensional square lattice 
is shown. 

The calculation for the clusters in Fig. was the 
largest possible. To have an additional comparison, we 
also performed calculations keeping only the operators 
which affect at most four different sites on the clusters of 
Fig. In summary, the various truncation schemes are 

(i) minimal model: only H(j, H t , £/nn 

(ii) NN: cluster used in 2nd order perturbation theory, 
Fig. Eli. 

(iii) plaquette and three-site chain: cluster used in 4th 
order, Fig. [5p. 

(iv) operators on up to four different sites on the clus- 
ters of Fig. 0i. 

(v) double-plaquette and four site chain: clusters 
used in 6th order, see Fig. 0:. 

The NN truncation (ii) yields the leading order result 
for the NN spin coupling Ji, the plaquette truncation 
(iii) yields the correct result in 4th order, and the double 
plaquette truncation (v) yields the correct result in 6th 
order. We confirmed that our implementation complies 
with these checks. 

The computation of the effective model splits into two 
parts. First, one has to set up the flow equation, i.e., 
one has to calculate the commutator [r],H]. According 
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Calculation 


^operators in H 


# of terms in DE 


(i) minimal 


4 


3 


(ii) NN 


6 


8 


(iii) plaquette 


172 


10,364 


(iv) 


2,217 


1,341,736 


(v) double-plaquette 


26,251 


304,514,721 



TABLE II: Number of operators in the Hamiltonian reduced 
by the symmetries given in the main text for various approx- 
imation schemes. This number coincides with the number of 
differential equations. In the last column the number of bi- 
linear terms on the right hand side the differential equations 
is listed. 



For the double plaquette, about 65 separate runs were 
performed, each with a CPU time below 10 hours and 
less than 0.5 Gbyte memory on Sun UltraSPARC work- 
stations. 

The second program solves the flow equation, i.e., the 
set of differential equations. As the set of equations is 
large this second program is very demanding both in time 
and in memory space. Typically 120 hours of CPU time 
and about 2.5 Gbyte are needed on Sun UltraSPARC 
workstations. 



IV. RESULTS 



to the truncation scheme newly generated terms are dis- 
carded or included in H, in rj, and in the calculation of 
[q, H]. This is done in a dynamic fashion in a CH — h pro- 
gram. The list of terms considered grows in the course of 
the automatized computation of the commutators. Since 
we consider operators on a finite number of sites there 
is a point, when no new terms are generated. At this 
point the set of differential equations, the flow equation, 
is closed. 

Second, one has to solve this flow equation. The dif- 
ferential equations have to be set up only once for each 
truncation scheme. Then they must be solved numeri- 
cally for each set of initial conditions given by the initial 
value of t/U. Since all coefficients are entangled in the 
differential equations the numerical integration provides 
all the coefficients - magnetic exchange couplings, hop- 
ping and interactions of charge - of the effective model 
together. 

To keep the number of terms as small as possible, we 
made use of the point-group symmetries of the square 
lattice, i.e., rotations about multiples of 7r/2 and reflec- 
tions about the axes and diagonals. Additionally, terms 
which are related by hermitian conjugation, by particle- 
hole transformation, and by flipping a —* —a all spin in- 
dices carry the same coefficient. Table UTI shows the num- 
ber of terms, reduced by the symmetries, in the Hamil- 
tonian of each calculation. As rj and H are both linear 
in the coefficients of the terms the right hand side of the 
differential equations is bilinear in these coefficients. The 
third column of Tab. ITT1 lists the number of such bilinear 
terms in the differential equation. The number of vari- 
ables and the number of differential equations is growing 
fast with the size of the cluster defining the truncation 
rule. This renders both the set up of the flow equation 
and its numerical integration computationally costly. 

We implemented the CUT calculation on the computer 
with two C++-programs. The first sets up the differen- 
tial equation by symbolic calculation of the commutator 
of rj and H. For a large number of terms in the hamil- 
tonian this task is very time-consuming. But it is possi- 
ble to split the calculation into a number of independent 
parts that can be accomplished by separate programs. 



Here we present the results for the five different trun- 
cation schemes (i) to (v) described above in Sect. 1111 Bl 
First, we discuss some properties of the continuous uni- 
tary transformation which enable us to analyze the limi- 
tations of the construction of the t-J model. Second, data 
on the effective magnetic couplings will be given. Third, 
data on the effective charge couplings will be shown. 

A. Properties of the CUT 




l-U 

FIG. 6: Evolution of the residual off-diagonality in the double- 
plaquette calculation (v) for various values of W/U. 



The CUT is induced by operators that change the num- 
ber of DOs. They are off-diagonal in the number of DOs. 
We define the residual off-diagonality as the sum per site 
of the squares of all coefficients that belong to operators 
that change the number of DOs. In the calculation the 
Hamiltonian is written in terms of products of the local 
operators in Table [i] This implies the definition of the 
coefficient as the prefactor of such a product. 

The residual off-diagonality measures the extend to 
which the unitary transformation has eliminated the 
terms that change the number of DOs. In a numeri- 
cal integration of the differential equation the residual 
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off-diagonality will not vanish exactly. Therefore, we de- 
fine a small but non-zero value and stop the flow when 
the residual off-diagonality falls below this value. Fig. 
displays the evolution of the residual off-diagonality 
in the double-plaquette calculation (v) during the flow, 
parametrized by the variable £U, see Sect. lIIII for various 
values of W/U. At I = the operators in T + 2 and T_ 2 
in J2J contribute to the residual off-diagonality. Written 
with the operators in Table |U these are 4*2 = 8 terms 
(hopping forth and back in x- and y-direction, with spin 
up or down) in T +2 and T_ 2 each. Therefore the curves 
in Fig. for the residual off-diagonality start at I — 
with the value 9>t 2 +2 + 8*-2 = 16 * 2 - 

For small values of W/U the residual off-diagonality 
shows exponential decay. For larger values of W/U the 
convergence is slower. For W/U = 1.12 it shows for the 
first time pronounced non-monotonic behavior. The dot- 
ted black curve, corresponding to W/U = 1.12, falls be- 
low 10 -27 to rise again up to 10 -3 for I ■ U w 270. Non- 
monotonic behavior occurs for all values W/U > 1.12; 
it can become a problem since we have to compare the 
residual off-diagonality to a small but finite value to de- 
cide at which value of I to stop the integration of the 
differential equations. But the values of the dominant 
coefficients of the effective model are reached already at 
small I. They are not affected noticeably by different 
cutoffs for the residual off-diagonality. 

Most importantly, one has to check whether the map- 
ping to the effective model is possible at all. Below we 
will calculate the apparent charge gap A app within the 
effective model. This is the minimal energy of a DO, e.g. 
a hole, propagating through a paramagnetic spin back- 
ground after the CUT. Since the paramagnetic spin back- 
ground is not the ground state the apparent charge gap 
A app is not the true charge gap. The apparent charge 
gap measures the separation of energy scales. This sep- 
aration governs the convergence of the mapping as we 
perform it. We will discuss this issue and its physical 
significance in more detail in the Discussion. 

We find that the apparent charge gap A app becomes 
negative (!) for values of W/U > 0.92. Once excitation 
energies become negative the CUT does not work because 
the limit I — * oo does not exist any longer. The #DO- 
changing processes can no longer be eliminated because 
they constitute no longer vacuum fluctuations around a 
stable vacuum. The negativity of excitation energies sig- 
nals that the chosen reference ceases to be a physically 
reasonable reference. 

The non-monotonic behavior of the residual off- 
diagonality appears only in a parameter range where 
the mapping is no longer possible. For parameters 
W/U < 0.9, where the mapping is possible according to 
A app > 0, we find an exponential decrease of the resid- 
ual off-diagonality down to values where the calculation 
is limited by the numerical accuracy of the computer im- 
plementation. A pronounced increase of the residual off- 
diagonality was only found in the double-plaquette calcu- 
lation (v). The calculations (i)-(iv) show an exponential 



drop of the residual off-diagonality for all values W/U . 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



W/U 

FIG. 7: Distribution of coefficients as function of W/U accord- 
ing to the spatial extension d of the respective terms in the 
double-plaquette calculation (v). The sum of the squares of 
the coefficients are plotted relative to the second order Heisen- 
berg exchange Jj 2 ' = 4i 2 /(7. 

Next, we analyze the distribution of coefficients over 
operators of various extensions d, see Sect. Mil and the 
effect of truncating the proliferating terms. Fig. shows 
the sum of the squares of all coefficients d of terms with 
the definite extension d at the end of the flow. They are 
given as functions of W/U. The results are given in units 
of j[ 2) = 4t 2 /U which is a natural scale for the parame- 
ters defining the effective model. The only operator with 
extension zero is the Hubbard repulsion term Hjj . For in- 
creasing W/U the relevance of the most local terms with 
d = and d = 1 decreases, whereas terms with larger 
extension become more and more important. Note es- 
pecially the increase of the coefficients of the operators 
with maximal extension d = 3. 

The distribution of coefficients indicates that up to 
W/U w 1.2 all important operators are included in our 
calculation. Beyond W/U > 1.2 the coefficients with 
extension d = 3 become more important. Thus terms 
with larger extension would have sizeable coefficients 
and should be taken into account. This indicates that 
the approximations become definitely insufficient beyond 
W/U > 1.2. 

To obtain clearer information on the applicability of 
the mapping and on the accuracy of the approximations 
we examine the apparent charge gap A app . To this end, 
the dispersion of a single DO moving in the para- 
magnetic spin background is determined by a Lanczos 
calculation 25 26 . In this calculation we treat the oper- 
ators as vectors as in the projection technique 27 . The 
effective Hamiltonian acts by commutation as a super- 
operator on the operators. This is conveniently denoted 
by the Liouville operator C := [H, ■]. The "scalar prod- 
uct" is defined by Tr(A^ Bpo) relative to the para- 
magnetic reference ensemble ©. Note that this scalar 
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FIG. 8: Apparent charge gap as calculated with Lanczos diag- 
onalization. Grey lines show the minimal eigen values, black 
lines the result of a linear extrapolation. 



product is only positive semi-definite since there can be 
operators whose norm vanishes. 

The projection formalism is the appropriate concept 
for this calculation since we are not dealing with simple 
states in a Hilbert space but with ensembles which are 
derived from the paramagnetic one © by the application 
of operators. These operators induce deviations from the 
reference ensemble which constitute excitations. As usual 
in the Heisenberg picture, the dynamics of operators in 
the frequency domain is captured by the commutation 
with the Hamiltonian. 

We start the Lanczos approach from the operator Oq 
which generates a single DO with momentum k 

where the number of sites is denoted by N and the sum 
extends over the whole lattice. A set of orthogonal oper- 
ators is generated by the iterative application of C 

O n+1 = CO n - a n O n - #0 n _i , (20) 

where 

On = {O n \CO n )/(O n \O n ) , (21a) 

b 2 n+1 = {O n+1 \O n+1 )/(O n \O n ), b = 0. (21b) 

Note that the application of C derived from the effective 
Hamiltonian after the CUT does not change the number 
of DOs. Hence its action on the initial vector is (i) to 
shift the DO and/or (ii) to change the spin background. 
The iteratively built basis {O n } describes a single DO 
(charge excitation) at momentum k including its mag- 
netic dressing. In this basis the Liouville operator is a 
tridiagonal matrix where the a n are the diagonal matrix 
elements and the b n are the elements on the secondary 
diagonal. 



The lowest energy in this subspace defines the (appar- 
ent) dispersion of a single charge excitation relative 
to the paramagnetic ensemble. Hence the evaluation of 
the lowest accessible energy in any truncated subspace 
of {O n } provides an upper bound to Eh- The apparent 
charge gap A app is finally given by 

A app = 2min£; fe , (22) 

k 

where the factor of 2 is put to account for the creation of 
a particle and a hole. For vanishing hopping A app = U 
holds. The apparent charge gap A app is not the true 
charge gap because the paramagnetic spin background is 
not the true ground state. But A app is a measure of the 
separation of energies between sectors of different double 
occupancy. 

Numerically, one has to restrict the above procedure to 
truncated basis sets. In practice, we construct a sequence 
of basis sets, labelled by the integer n, by applying the 
Liouville operator Co := [To, •] (cf. Eg . I4b|l iteratively to 
the initial operator (|19fl . The first basis n = 1 consists 
only of the single operator Oq. The subsequent basis 
sets are generated by applying T (n — 1) times. All the 
components which are products of the local operators in 
Tab. [I] are treated as independent vectors of the basis. 
For each basis set the lowest eigen value of the matrix 
representation of the Hamiltonian in this basis yields an 
upper bound to the lowest energy E). of a single DO. 

The calculation was done for k- values on the high sym- 
metry lines in k-space from (0, 0) — (ir/a, n / a) — {ir / a, 0) — 
(0,0). Figure |S| shows the results of the Lanczos calcu- 
lation. The computational effort grows with the number 
of terms in the Hamiltonian. For the effective Hamilto- 
nians (iv) and (v) it was only possible to do the calcu- 
lation up to the basis set n = 6, whereas for the NN 
(ii) and plaquette (iii) Hamiltonians the n = 7 calcula- 
tion was feasible. In Fig. [5] grey lines show the smallest 
eigen value found for all momenta k in the largest possi- 
ble Lanczos calculation. An additional estimate for the 
gap is found by extrapolating the lowest eigen values at 
momentum (jr/a,Tt/a) for each basis set as function of 
1/n to n — > oo. 

The lowest eigen value at momentum [it / a, it / a) is 
usually also the minimum of the dispersion except for 
n = 5, 6 or 7 where the whole dispersion becomes very 
flat. The result from a fit linear in 1/n is included in 
Fig. [S] by black lines. For the minimal eigen values as 
well as for the extrapolated ones the result for the Hamil- 
tonian (iv) is covered by the double-plaquette result (v), 
since the coefficients of the Hamiltonians nearly coincide. 
The lowest eigen values for the Hamiltonians (iv) and (v) 
are larger than the lowest eigen values for the Hamiltoni- 
ans (ii) and (iii) because larger values of n were accessible 
for the latter. 

A first guess for the gap is A app /[/ = 1 — W/U, see 
e.g. Ref. H- So one expects that the gap closes around 
W/U = 1. The apparent charge gap calculated within 
the effective model using Lanczos-diagonalization dis- 
plays a similar behavior. The gap decreases almost lin- 
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early from A app /U = 1 at W/U — 0. There is, however, 
a certain downward curvature so that it closes at about 
W/U ~ 0.9. Taking into account even larger basis sets 
might further lower the ratio W/U where the gap closes. 
The linear extrapolation corroborates the results from 
the linear eigen values to a certain extent. But we could 
not identify a clear asymptotic behaviour for large val- 
ues of n. Other extrapolation schemes, for instance as 
function of 1/y/n, point towards an earlier closing of the 
apparent charge gap. Therefore the value W/U » 0.9 for 
the closing of the gap must be seen as a rough estimate. 
The determination of the power law by which the gap 
closes is beyond the scope of the present paper. 

The above findings are discussed in more detail in Sect. 
Ivl In the sequel, we have to keep in mind that the map- 
ping to the effective model conserving the number of DOs 
is not possible for W/U > 0.9. For completeness, we will 
present the coefficients of the generalized t-J model up 
to values W/U = 1.6 to illustrate in which way the re- 
sults are affected by the breakdown of the mapping. As 
a reminder of the breakdown, the definitely unphysical 
region beyond W/U w 1 will be shaded in the figures. 



B. Effective spin model 

Strictly at half-filling, the effective model can be cast 
into the form of a pure spin Hamiltonian. There are 
various two-spin couplings of the Heisenberg-type 



^2-spin — J]] J I 

i,3 



\j- 



I Si ' S*j 



(23) 



This type of exchange may exist for sites i and j of any 
distance but it is of course largest for adjacent sites. 
Fig. displays the coefficients of various two-spin cou- 
plings. The Heisenberg NN exchange J\ is shown in 
Fig-UK; the exchange couplings between second and third 
nearest neighbors, J2 and J3, are depicted in Fig. [[J) and 
EJ, respectively. The coefficients are shown relative to 
j[ 2 ' = At 2 /U. The coupling J\ is only renormalized 

(2) 

slightly to lower values w 0.92 J) ' for W/U — 1 con- 
trary to what one might have expected from the strong 
influence of higher terms in perturbation theory. Already 
the NN calculation reproduces the leading order correctly 
and yields a very reasonable result for larger values of 
W/U, see also Fig. 0] and the discussion thereof. The 
convergence of the various calculations shows that the 
terms retained are sufficient to determine J\ to very good 
accuracy. Note that the dotted curve of the plaquette 
calculation (iii) is almost hidden by the double-plaquette 
result (v). 

The couplings J2 and J3 behave like {t/U) A for small 
W/U. Their absolute convergence is as good as for J\. 

But due to their small values - both are smaller than 

(2) 

0.02J{ for W/U < 1 - the relative differences between 
calculations on different clusters are still discernible. 

There are not only exchange terms between two spins 
but also between four or more spins. Exchange between 
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FIG. 9: Couplings of two-spin operators relative to J-J as 
function of W/U. 



an odd number of spins is not generated since such terms 
would change sign under time reversal. The first term 
that one encounters in perturbation theory 3 in 4th order 
is the four-spin part of the ring-exchange which takes the 
form 

H n = Ju Yl [(Si -Sj-XSfc- S { ) + (8^80(8^8*) 



(Si-SfcXSj-s, 



(24) 



where (i,j,k,l) label the sites around a plaquette of the 
square lattice in cyclic order. In 6th order an additional 
independent four-spin operator reads 



E 

(i,j,k,l) 



(Si-sots,- -S, 



(25) 



The coefficients J a and J x are plotted in Figs. 110b and 
110b . The ring exchange Jn is the most important mod- 
ification of the pure NN Heisenberg model. Its value 
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FIG. 10: Couplings of four-spin operators relative to J{ as 
function of W/U. 



(2) 

is larger than 0.2J] ; for W/U = 1. The relevance of 
the ring-exchange was also confirmed experimentally in 
the CuCvplanes of I^CuO^ and theoretically in the 
derivation of the effective spin model from a three-band 
model22i2L2I. The influence of ring-exchange on the Ra- 
man line shape was computed in spin wave theory^. 

It is worth noting that the coefficients of the effective 
spin model are smooth functions of W/U up to W/U = 
1.6. The vanishing of the apparent gap A app and the slow 
drop of the residual off-diagonality for W/U > 1 do not 
give rise to an anomalous behavior of the effective spin 
coefficients. The convergence of the results for values 
W/U < 1 is very good, i.e., the couplings do not change 
much as function of the cluster size for the larger clusters. 



C. Effective charge model 

The unitary transformation deals also with the motion 
and interaction of charges. In the course of the trans- 
formation charge excitations (DOs) are created virtually. 
After the transformation, DOs can be introduced exter- 
nally by doping or thermal fluctuations. Note that gener- 
ically the effective model displays a finite U term. So DOs 
are not suppressed. But due to the disentanglement of 
the sectors of different number of DOs there is no direct 
influence of these sectors on one another. For instance, 
the magnetism strictly at half-filling is not influenced by 
the physics of the sector with a particle-hole pair (two 
DOs in our counting). 

But thermal expectation values are generically influ- 
enced by the presence of sectors of different number of 



DOs. Also spectral properties are influenced by the pres- 
ence of these sectors. The spectral weight will be dis- 
tributed over a number of sectors. There are Hubbard 
bands beyond the lower (one hole) and the upper (one 
electron) Hubbard band, for instance a trans-upper Hub- 
bard band characterized by two electrons and a hole rel- 
ative to the paramagnetic reference ensemble. But pre- 
liminary estimates of the spectral weight in this band 
indicate that it is below the order of one percent. Yet it 
is a strong point of the approach chosen that the concept 
allows to discuss such subtle effects. 

In the present paper, we focus on the effective Hamil- 
tonian leaving the determination of spectral weights to 
future work. We will present the results for the most 
important terms concerning the dynamics of charges. In 
order not to be lost in a proliferating number of terms we 
will present results for terms which are at most of order 
t 2 /U in the limit of weak hopping. 

First, we address the terms which move a single DO, 
i.e. various hopping processes. Besides the NN hop- 
ping To, the next relevant hopping processes are those 
between next-nearest neighbor (NNN) sites and third- 
nearest neighbor sites (3NN). The NNN sites lie on di- 
agonally opposite corners of the plaquettes of the two- 
dimensional square lattice. All hopping processes and 
interactions that appear in second order in t/U will be 
shown. The NNN hopping Tq reads 

T o = l ' E [( X ~ n i,<r) c \,w c j,A l - n j.a) 



n i,A,a c j,a n j,a + h ' C - 



(26) 



where stands for NNN sites and a = —a. The 

spin-dependent NNN hopping T' s reads 



+ h/L^/jJ • S fc + h - c J (27) 



where (i, k,j) stands for three sites wherein i and j are 
NNN and A: is a nearest neighbor to both i and j. The 
symbol er = (a Xl a y , a z ) stands for the vector made out 
of the three Pauli-matrices; represents the usual S — 
1/2 spin vector at site fc. Hopping and spin-dependent 
hopping between 3NN is also produced in second order 
in t/U. They read 



(28) 
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where stands for 3NN sites and 

((i,kj));a,l3 1 

where ({i,k,j)) stands for three sites wherein i and j 
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FIG. 11: Effective hopping coefficients between (a) NN sites, 
(b) NNN sites, (c) spin-dependent NNN hopping, (d) 3NN 
sites and (e) spin-dependent 3NN hopping. 

are 3NN and k is a nearest neighbor to both i and j, so 
k is just the site between i and j. All these hoppings 



(|26[1 - <|29[) conserve the number of DOs. Figs. Illb-e show 
the coefficients of the NN hopping to (Tq from Eq. I |4bjl ). 
NNN hopping t', spin-dependent NNN hopping t' s , 3NN 
hopping t" and spin-dependent 3NN hopping respec- 
tively. They are shown in units of the unrenormalized 
value to(0) = t for the NN hopping. 

The NN hopping to shown in Fig. Illb remains un- 
changed from its initial value for the NN truncation. In- 
cluding more terms reduces the effective to- At W/U ~ 
1.36, the double plaquette calculation (v) shows a sudden 
change of slope. This is related to the non-monotonic 
behavior of the residual off-diagonality. But the value 
W/U fa 1.36 is already beyond the point W/U w 0.9 
where the apparent gap A app becomes negative and the 
mapping is no longer possible. The signs of t' (NNN hop- 
ping) and of t" (3NN hopping) are opposite to the sign 
of t (NN hopping) . The value of t' is approximately two 
times t" because it is generated from the two possible 
hopping processes to go from one corner of a plaquette 
to the opposite corner. In contrast, there is only one pro- 
cess to generate t" in second order. The results for the 
calculations (iii)-(v) lie almost on top of each other for 
the coefficients t, t' and t" until spurious behavior sets in 
at larger values of W/U. 

The coefficients of the spin-dependent hoppings t' s and 
t" are shown in Fig. II lb and e. They increase quadrati- 
cally for small W/U just like t' which leads to the linear 
behavior in units of to depicted in Fig. These hop- 
pings are terms which concern three sites. Their coeffi- 
cients have approximately the same value since it does 
not matter very much whether the three sites arc aligned 
linearly or form an angle of 90 degrees. We emphasize 
that both coefficients t' s and t" are of similar size as t'. 
This implies that the commonly used phenomenological 
description in terms of a t-t'-J model is not justified. If an 
extension of the simplemost t-J model is used the NNN 
and 3NN hoppings, t' and t", and the spin-dependent 
hoppings t' s , t" s should all be included. 

Figures IT21 and H31 show the coefficients for terms that 
describe interactions between charges. In Fig. I12h . the 
Hubbard repulsion U is shown in units of its unrenormal- 
ized value. All calculations show a slight increase of the 
effective U with W/U. It is interesting to note that this 
basic quantity neither diverges nor stays constant in the 
process of the disentanglement of the various sectors of 
double occupancy. 

All processes discussed in the sequel have the property 
to be active only if at least two DOs are present. In our 
counting this means that at least two holes or two doubly 
occupied sites or one hole and one doubly occupied site 
must be present. In this sense, the following processes 
represent two-particle interactions. The density-density 
interactions for NN sites is described by 



(30) 



where n* 



n i,1 + n i,l 



1. Similarly, the operator for 
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FIG. 12: Hubbard repulsion U (a) and NN interaction V (b) 
and NNN interaction V' (c) between holes or double occu- 
pancies. 



density-density interactions for NNN sites reads 



V n ^ n r 



(31) 



The coefficients V and V' are shown in Fig. 112b and 
112b . They are given in units of the second order scale 
j[ 2 ^ = At 2 /U . The attractive interaction V is propor- 
tional to — t 2 /U = —J^/4 as expected from perturba- 
tion theory. It represents a sizeable contribution at all 
W/U. In contrast, the V coupling NNN sites, which is 
of fourth order in t/U, turns out to be very small 33 . This 
fact justifies our choice to focus otherwise on the terms 
which appear already in second order in t. 

There are two additional types of operators that are 
created in second order in t/U. The first type consists 
of operators which destroy or create two electrons at the 



same site. Figure ITSk shows the coefficient Vp a ir of 



-^pair — ^pair 



i c j,i c j,r 



(32) 



Figure IT3b and c show the coefficients V^ air and V^' air of 

^p'air E [Ws^Vj,^ 1 ~ n 3,^ 
(i,k,j),cr 

+4 )0 4,a C i,<r( 1 - n i,JS>^ + h - C - ' ( 33 ) 



pair 



H, 



pair 



^pair E [4,<r C k,a C i,a n i,<r C jA 1 ~ n j,^ 



+ C L C L C i,a( 1 



r)c,<x n 7.<T +ll - C - 



(34) 



The operators -ff pair and i?p' air describe processes where 
two electrons from sites i and j are transferred to site 
k. They do not change the number of DOs since at one 
site (with local operator Ci. a (l — n^g-)) a DO is created 
while at another [ci^n^g) a DO is annihilated. The local 
operator c^c^ , turns an empty site into a doubly occu- 
pied site and thus does not change the number of DOs in 
our counting. Both the coefficients V' w and V'' ir change 
only slightly on passing from calculation (iii) to calcula- 
tion (iv); the difference between (iv) and (v) is minute 
so that we consider the final result as reliable for not 
too large W/U where the mapping to the generalized t-J 
models is possible. 

The last type of operators created in second order in 
t/U are correlated hopping terms. We classify them as in- 
teractions because they have a non-vanishing effect only 
if at least two double occupancies or holes are present. 
The operator 



v' = v; 



E 

(i,k,j);a,f3 



[» 



(1 



- n i, a 4,ci c jt n j,0 n k + h - c - 



(35) 



describes a process of a DO hopping between the NNN 
sites i and j if and only if there is a DO on the site k. 
The corresponding operator for a 3NN process is 



v n,0 v n 



E 

((i,k,j));a,p 



(! - ^i,a) C U C ^( 1 " n j,p) n k 



' n i,a c i,a C j^ n j,/3 n k 



h.C 



(36) 



The coefficients V/ t and V// are shown in Fig. I13H and e. 

The above results define the charge dynamics of the 
t-J model in a systematic and exhaustive way. We high- 
light the importance of spin-dependent hoppings. All the 
values for the hopping coefficients and the interactions 
show excellent convergence in the relevant region up to 
W/U w 0.9. We find again that for values W/U > 0.9 
the coefficients display anomalous behavior. This ef- 
fect is most striking for the coefficient V' in Fig. 112b . 
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We attribute this spurious behavior to the breakdown of 
the mapping of the Hubbard model to a generalized t-J 
model without specifying the state of the spin degrees of 
freedom. 



V. DISCUSSION 

In this section, we will discuss the physical significance 
of our findings for the square lattice and compare them 
to the situation for the Hubbard model in infinite dimen- 



In the derivation of the effective t-J-model, one does 
not aim at the solution of the complete problem of the 
Hubbard model on the square lattice. The CUT is de- 
signed to disentangle the energetically low-lying dynam- 
ics, i.e., the dynamics of the spins and of the doped holes, 
from processes involving more double occupancies. So no 
particular state in spin space is assumed, but an effective 
model without processes changing the number of double 
occupancies is derived. In this calculation, the paramag- 
netic ensemble at half-filling serves as reference ensemble. 

The general consideration presented in the Introduc- 
tion suggested that the reduction to the effective model 
without charge fluctuations cannot be defined for too 
large W/U. Our findings fully corroborate this view. 
The reduction of the Hubbard model to a generalized 
t-J model is possible as long as the apparent gap A app 
remains non-negative. The appearance of negative values 
for A app is an artefact. It indicates that the paramag- 
netic reference ensemble does no longer represent a phase 
which is stable against fluctuations involving double oc- 
cupancies. In this case the reduction to an effective model 
without virtual double occupancies is no longer possible. 

The concept of the apparent gap is introduced as a 
quantitative measure of the energetic separation of the 
low-lying degrees of freedom from the fluctuations involv- 
ing double occupancies. But it is not easy to visualize 
since it does not refer to true eigen states of the system. 
So it is helpful for a qualitative understanding to consider 
a modified system where the apparent and the true gap 
coincide. This situation is realized for a Hubbard model 
on the Bethe lattice in the li mit of infinite coordination 
number Z — > oo, cf. Refs. For infinite dimen- 

sional lattices it is assumed that long-range order can be 
completely suppressed by frustration, e.g., in the gener- 
alizations of fee lattices to infinite dimensions'^. Since 
the magnetic couplings J scale with the inverse coordina- 
tion number Z~ x there are no magnetic correlations at all 
once the static sublattice magnetization is suppressed 37 . 
This is so since the second moment of the magnetic cou- 
plings, averaged over adjacent sites, vanishes in the limit 
Z — > oo. There is even no nearest- neighbor correla- 
tion. Hence the paramagnetic reference ensemble rep- 
resents the highly degenerate magnetic ground state in 
infinite dimensions without long-range order. The appar- 
ent gap and the true charge gap are identical. Their clo- 
sure signals a real insulator-to-metal transition. In finite 
dimensions the closure of the apparent gap signals only 
a crossover from well-separated to non-separated energy 
scales. 

We see that the (paramagnetic) insulator-to-metal 
crossover in finite dimensions can be studied in the puri- 
fied form of a phase transition in the infinite dimensional 
Hubbard model. Indeed, similar behavior is found. The 
single-particle gap A vanishes in infinite dimensions at 
about 1.1W - 1.2VF 10 i 34 i 38 i 39 i 40 . It disappears roughly 
linearly as function of LSi2i2£. These similarities corrobo- 
rate the infinite dimensional model as illustration for the 
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crossover in finite dimensions. 

The validity of our approach is clear for small values 
of W/U. The only approximation used is the truncation 
of processes of a spatial range d > 4. The remaining 
question is to know up to which ratio of W/U the results 
are reliable. Our findings in finite dimensions illustrate 
that the answer to this question is difficult already on the 
conceptual level. The apparent gap A app measuring the 
separation of energy scales does not compare the energies 
of real eigen states. Hence we claim that the validity of 
the mapping of the Hubbard model to a generalized t-J 
model is not limited by a sharply defined transition but 
by a gradual crossover. 

An estimate up to which ratio of W/U the mapping 
is reasonable can be obtained from the closure of the 
apparent gap at around W/U ~ 0.9. A better estimate 
for the region of validity should compare A app to the 
size of the coefficients of the neglected operators. But 
these coefficients are not available. Hence we compare 
A app to the size of the coefficients of the operators with 
the maximal extension considered, that is d — 3. Both 
energies become equal for W/U « 0.85. If in addition 
we take the uncertainty of the extrapolation procedure 
into account we arrive at the conservative estimate that 
the generalized t-J model can be used up to W/U » 0.8 
which includes the commonly assumed parameters for 
cuprate planes t/J » 3 which translates to W/U m 0.7. 
Up to W/U « 0.8, the quantitative results provided in 
this article are reliable. 

In conclusion, we find that the possibility to map 
a Hubbard model in a systematic and controlled way 
to a generalized t-J model depends essentially on local 
physics. We stress, however, that the breakdown of this 
mapping does not represent a real physical phase transi- 
tion in finite dimensions. We suppose that the breakdown 
of the mapping manifests itself in physical properties as 
a crossover. It is plausible, for instance, that the nature 
of the elementary excitations changes: holes and mag- 
netic modes for large values of the interaction become 
fermionic quasi-particles at small values of the interac- 
tion. 



VI. SUMMARY 

We introduced a self-similar CUT to map the Hub- 
bard model at strong coupling onto an effective model. 
The generator of the transformation is chosen such that 
the number of double occupancies is conserved by the 
effective model which is thus a generalized t-J model. 
The CUT allows to do this in a systematic and non- 
perturbative way which includes the dynamics of doped 
holes as well as the spin dynamics. Thus our analysis ex- 
tends previous ones in two important aspects: the non- 
perturbative treatment of the couplings and the system- 
atic treatment of the dynamics of doped charge carriers. 

The calculations confirmed our expectations which 
were based on qualitative arguments. The fundamental 



result is that for U > 1.3W the mapping of the Hubbard 
model to a generalized t-J model is possible. The param- 
eters of the effective model can be determined reliably by 
the CUT. They are governed mainly by local physics. On 
the other hand, the reduction of the Hubbard model to 
an effective model without charge fluctuations is not pos- 
sible for U < 1.2W. Here, the sectors of different number 
of double occupancies cannot be disentangled. 

The truncation scheme necessary to close the flow 
equation of the CUT was based on the extension of the 
operators. In order to define the extension of an operator 
a physically meaningful unique notation is required. This 
was achieved by introducing a normal-ordering with re- 
spect to a vacuum. The vacuum chosen here was the 
paramagnetic ensemble without any magnetic correla- 
tions. 

Besides the NN Heisenberg exchange further two-spin 
terms and four-spin terms were taken into account. It 
was shown that two-spin couplings other than the Heisen- 
berg term are small whereas the four-spin ring exchange 
Jo gives a sizeable contribution to the effective spin 
model. Results for the size of the coefficients for NN hop- 
ping, NNN hopping and spin-dependent NNN hopping 
were obtained. It was found that the spin-dependent 
NNN and 3NN hopping is as important as the spin- 
independent one. We recommend to account for this 
fact in phenomcnological parametrizations of experimen- 
tal Fermi surfaces. In addition, the size of the interac- 
tions between holes and double occupancies on NN and 
NNN sites were calculated. Various truncation schemes 
show very good convergence of the coefficients of the ef- 
fective model in the relevant parameter region W < 0.8J7. 

Our results suggest that self-similar continuous unitary 
transformations can be used also for other models, e.g. 
more realistic multi-band Hubbard models, in order to 
achieve a systematic and quantitative reduction to effec- 
tive models. 
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APPENDIX A: THE NN TRUNCATION 

To carry out the CUT with all operators on NN sites 
the terms l|12c[l and <|12d(l in equation (|12fl have to be 
included in the Hamiltonian 

H v = V{l)^nifij (Al) 

tfpair = V paiT (£) 4Al C U C i,1 ■ ( A2 ) 

(i,3) 

Hy describes interactions of electrons on NN sites, i? pa ir 
the hopping of two electrons between site i and j. 
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These terms do not change the number of DOs. Thus, 
they generate no new contributions in r/. Calculating 
[q, Hnn + Hy + //pair] an d neglecting again terms that 
do not fit on NN sites we obtain the closed flow equation 



d e U(£) = 32t +2 (£) 2 (A3a) 

d e t (£) = (A3b) 

d t t +2 (l) = -2U(£)t +2 (£) (A3c) 

+t +2 (£) (2V{£) - 2V pair {£) - 3J(€)/2)) 

d t Ji{t) = m +2 (£) 2 (A3d) 

d e V{£) = -U +2 {£) 2 (A3e) 

^VpairW - 8t +2 (£) 2 . (A3f) 



The solution for J\ is 



where r = y28tg + £/ 2 and D = arctanh([/o/r). This 
yields in the limit £ — > oo the effective NN spin-coupling 
of the NN truncation 

Jf&^VlF + MP-jU, (A5) 

where the final result is given in the bare parameters 
recalling to = t, Uo = U. 
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